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Abstract 

A general study of excitations up to two-phonon states is carried out using 
the intrinsic-state formahsm of the Interacting Boson Model (IBM). Spec- 
tra and transitions for the different dynamical symmetries are analyzed and 
the correspondence with states in the laboratory frame is established. The 
influence of multi-phonon states is discussed. The approach is useful in prob- 
lems where the complexity of the IBM spectrum renders the analysis in the 
laboratory frame difficult. 
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I. INTRODUCTION 



Whether double-7 excitations exist or not in nuclei has been under discussion during the 
last 30 years. Double-7 excitations correspond to K'^ = 0^ and K'" = 4+ states, and are 
the bandheads of two rotational bands superimposed on the vibrational spectrum. Recent 
experimental improvements in nuclear spectroscopy following Coulomb excitation |Q, in- 
elastic neutron scattering and thermal-neutron capture have made possible the study 
of highly excited low-spin states which until now had been inaccessible. In transitional 
even-even nuclei i86-i92Qg ^^^pt large components of the double-7 excitation have been 
identified In well-deformed nuclei the search for double-7 excitations has led to the 

identification of a first candidate namely the K"^ = 4"*" level at 2055 keV in ^^'^Er Sub- 
sequently, other two-phonon excitations have been observed. For example, based on their 
decay properties the levels in ^^^Er at excitation energies of 1943 and 2028 keV were judged 
to be collective and interpreted as = and K'^ = 4+ members of the double-7 mul- 
tiplet, respectively 0J^. Other examples of double-7 excitations are the K'^ = 4+ state at 
1435 keV in ^^^^Mo |], K"" = 4+ states in ^^^-^^^Gd |10| or if^ = 4+ state at 2173 keV in 
^^"^Dy ||lT|. It should be noted, however, that the identification of such states as members 
of a two-phonon multiplet is still very much under debate ||12|| . 

A possible framework for investigating two-phonon excitations is the Interacting Boson 
Model (IBM) |]13|. This model has been very successful in describing low- lying collective 
spectra in medium-mass and heavy nuclei. Two-phonon excitations are present in the IBM 
but the analysis of their properties is not so easy to carry out in the laboratory frame since 
in the energy region of interest (around 2 MeV) the density of states is high and their decay 
properties are intricate. For example, it is as yet not clear how and, if so, to what extent 
anharmonicities can be introduced in the model, a question that did generate controversy 



in the early days of the model [14,15 



In this paper it is shown that the intrinsic-state formalism p6| , [T7[] provides an ideal frame- 
work for analysing the problem of two-phonon excitations in the IBM. The intrinsic-state 



2 



formalism is an approximation (of order 1/N) to the exact IBM, but its advantage is that it 
gives a clear interpretation of the structure of the states. In order to address the problem of 
two-phonon excitations we present here calculations in the Tamm-Dancoff Approximation 
(TDA), extended up to two-particle-two- hole excitations carried in the intrinsic frame of 
the IBM. This approximation is useful due to its simplicity and its physical transparency, 
and it allows to derive analytical expressions for energies and electromagnetic transitions 
applicable to a wide range of situations. 

The structure of the paper is as follows. First, a brief recapitulation of the intrinsic- 
state formalism as applied to the IBM is given, which is subsequently extended to include 
up to two-phonon excitations. Problems related to spurious states and anharmonicity are 



discussed. In section |TT| the results of schematic calculations in the intrinsic frame are 
compared with those in the laboratory frame for the three dynamical symmetry limits of 
the IBM. The character of the two-phonon states is established. Section |^ is devoted to 
electromagnetic transitions. Finally, in section [V] conclusions of this work are presented. 



II. THE IBM IN THE INTRINSIC FRAME 

The IBM describes low-lying collective excitations in even-even nuclei in terms of bosons 
with angular momentum (s bosons) and 2 (d bosons) |jl3|]. These bosons interact via a 
hamiltonian that is rotational invariant and number conserving and usually includes up to 
two-body interactions, although higher-order terms have been sometimes included. The 
most general two-body IBM hamiltonian may be written in multipole form as 

4 

H = Sshs + Sdhd + KoP^ P+T. ^lTl ■ Tl, (1) 

L=l 

where 

pt = l^t . _ l^t . (2) 
2 2 ^ ^ 

and 
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Tlm = Yl xie^hl X 7^2 )lm, (3) 
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The symbol ■ denotes scalar product, in this work the scalar product is defined as Tl ■ Tl = 
J2{~^)^ '^'^TlmTl-a4, while 7 with i = {£ = 2) corresponds to an s boson {d boson). 

M 

In Eq. d^) the operator = (— l)^~™'7^-m is introduced so as to verify the appropriate 
properties with respect to spatial rotations. 

The intrinsic-state formalism []TB|,|T7| provides a connection between the IBM and the 



Bohr-Mottelson collective model Let us consider a system consisting of a large but 

finite number N of interacting bosons. The dynamical behaviour of this system can be 
described in lowest-order as arising from independent bosons moving in an average field. 
In this approximation, the ground state of such a system is a condensate |c) of bosons, all 
occupying a single state Fj of lowest energy, 

\c) = ;^(rD^|o), (5) 



where 



= TrTF {'^ + ^ ^0 + sin 7 (4 + dU)^ . 



(6) 



The dynamical variables (3 and 7 are related to the quadrupole shape variables of the ge- 
ometrical model [0 (/3 > 0). The expectation value of the IBM hamiltonian (|I]) in the 
ground state (H) is 



c\H\c) = — — f 5 £s + 25 ^2 + (5 - 3 Ki + 5 K2 + 5 «:2 - 7 + 9 K 

^ (35 Ko + /?2(-70 Ko + 560 K2) - 80 VU 0^ x cos(3 7) K2 



5(1 + 



140(1 + /32 

+/?^(35/€o + 40x'fi:2 + 72/€4))- (7) 



Alternative expressions can be found in Refs. |16tp^p0[ using normal-ordered liamiltonians. 



The equilibrium values of the deformation parameters are obtained by minimizing the energy 
expression (|^ with respect to (3 and 7. Due to physical considerations a minimum for (3 00 
must be excluded I^Oj. Furthermore it follows from that a triaxial minimum is excluded 



and only prolate (7 = 0°) or oblate (7 = 60°) solutions are allowed. For an attractive 
quadrupole force < 0) these are found for x < and x > 0, respectively. A special 
situation of 7-unstability occurs if the quadrupole term vanishes (/t2 = 0) or if x = 0. 

The condensed boson V\ (P) is the first component of a new boson basis (deformed 
or intrinsic bosons). These deformed bosons are related to the spherical ones (laboratory 
bosons) through a general unitary transformation 77, 

£m P 

the deformation parameters r^f^ verifying the orthonormalization conditions 

^V^ilvlm = ^PP'i J2Ve*mV£'m' = ^tl'^mm' ■ (9) 

P 

The index p labels the different deformed bosons. After obtaining the r/ parameters for 
•p = cii.e. the equilibrium deformation parameters), the r] parameters for the excited bosons 
follow from the orthogonality conditions (^). An appropriate choice of deformed bosons can 
be found in Ref. |2D|: 

1 ^ /3cos7 ^ ^ /5sin7 
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7+ • 7+ 7+ COS 7 

mo = - sm 7 , r/22 = ^2-2 = 





1 




1 




~V2 ' 


V2-1 = 


V2 


V 

V21 


1 

~V2 ' 


V 

V2-1 = 


-1 
V2 


7— 
'?22 


1 

~V2 


7~ 
^2-2 — 


-1 
V2 



(10) 



other coefficients being zero. This basis has a clear correspondence with the usual (3 and 7 
excitations. 



In the following section the usefulness of this notation for the calculation of matrix 
elements in a one- and two-phonon basis shall become clear. 



A. One- and two-phonon states 

The excited bands can be considered as intrinsic excitations built on the boson condensate 
(I), obtained by replacing a c boson by an excited one, two c bosons by two excited ones 
and so on. This set of excitations contains spurious states that are coupled to the true 
physical ones and must be removed [^. This problem is a consequence of TDA. In the 
Random Phase Approximation (RPA) all the spurious states come out with zero energy and 



are decoupled from the physical states |2T|. However, RPA is more difficult to implement 
than TDA and obscures any physical interpretation. The problem of the spurious states will 
be returned to later. 

The excited states included in the present analysis are 

|p) = , ^ rt(rt,)^-'|o). 

, (11) 

where p, 7^ c. Eigenf unctions are linear combinations of these excited states plus the 
ground state, 

= iV« lw^\c) + E^^pb) + T.y^PP'\PP')] , (12) 
\ p pp' J 

where is a normalization constant. The ground-state contribution comes from the off- 
diagonal matrix elements of the hamiltonian between the ground state and the two-phonon 
states. 

To obtain the states (p!^), the hamiltonian is diagonalized in the basis of excited states 
(pup plus the ground state (^). The different matrix elements that intervene are 

{c\H\p) = VnF^I:> + 2v^(iV - l)Fi4 = 0, (13) 
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mw) = ^^/^' fZ, (14) 



{p\H\p') = F^;,) + 4(iV - + 5ppiN - 1) (f« + (N - 2)Fi£) , (15) 



{P\H\P'P") = I , , (4i^Sy' + V^ip'^' + 

+ 2(iV - 2)(5,,,Ff)^„ + V^c?i')) , (16) 



1 + 6pipi' 



(^pp'\H\py") = V' + ^P-W _ 2)^1) + (AT - 2)(iV - 3)i^£) 

i + OppiOppiiOpiipii, ^ ' 



+ dpp'){l + Vp'") 
+ V" (^p'p" + 4(iV - 2)F^X'c) + Vp" (^i?" + 4(iV - 2)Ff ) 



cpp"'c 



+ W' {Kv" + - 2)Fi;)„,) } , (17) 



where 



PpJ - XI ^^l^fimi^fimi (18) 



pp 

iimi 



and 



F(2) _ T/ „Pl* „P2* P3 P4 /-|Q\ 

P1P2P3P4 ''^^imi/2m2/3m3/4m4 7£imi7^2m2 '^3^3 '^4ni4" V^/ 

^imi^2™2^3™3^4'Tl4 

The coefficients include the single-particle energies ee plus contributions from the two- 
body terms in the IBM hamiltonian (|l]). The coefficients V£jmi/2m2/3m3/4m4 are defined 
as 

^imi,£2"l2/3"l3/4"l4 ~ mim2'\/l ~^ ^liln^mo,mi • 

(20) 

By construction F^, is symmetric under interchange of p and p' and -^pip2P3P4 same 
symmetry properties as Vfimi.^ama ,€3m3Am4- 
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B. Spurious modes. Goldstone bosons 



Before diagonalization spurious bosons must be removed from the basis. The procedure 
to separate out those modes from the physical ones is described in detail in Ref. [^. In 
this section we just enumerate the spurious (Goldstone) modes that occur in the different 
geometrical limits of the IBM. 

• P = 0, spherical nuclei. 

— 0(3) or 0(5) scalar hamiltonian. There are no Goldstone bosons. 

— 0(6) scalar hamiltonian. A spherical nucleus corresponds to the U{5) dynamical 
symmetry, not compatible with a symmetry group 0(6). Therefore this possibility 
must be excluded. 

• /3 > 0, deformed nuclei. 

— 0(3) scalar hamiltonian. For 7 = 0, F^. and F^ are spurious bosons. For 7 = 7r/3 
the spurious bosons are F^. and F^_. In the triaxial case F^., F^^ and F^_ are 
Goldstone bosons. 

— 0(5) scalar hamiltonian. All the excited bosons are spurious except the F^ boson. 

— 0(6) scalar hamiltonian. All the excited bosons correspond to Goldstone bosons. 



C. The anharmonicity problem 

It is not clear a priori to what extent f3 and 7 vibrations are harmonic in the IBM. The 
analysis of this problem is quite difficult in the laboratory frame; with the intrinsic-state 
formalism, however, some general conclusions can be obtained for arbitrary hamiltonians. 

Let us consider a system with a large number of bosons A^. For large A^ the IBM 
hamiltonian is approximately diagonal in the basis of excited states (pT]). The degree of 
anharmonicity q can be defined as (with this definition the harmonic limit corresponds to 
^ = 0) 
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with 



E;^ = {p\H\p)-{c\H\c), (22) 

E;^ = {pp\H\pp) - {c\H\c) . (23) 

These can be evaluated with the help of the following expressions, in which a Hamiltonian 
including up to three-body interactions have been supposed, 

{c\H\c) = NiN - l)(iV - 2)f£L + NiN - l)Fii + NF^\ (24) 

{p\H\p) = (AT - l)(iV - 2)(iV - 3)Fi£,, + (iV - l)(iV - 2)F(£ + 9(iV - 2)(iV - 3)Fig^,, 

+ (iV-l)(F«+4F(,y (25) 



(pplffbp) = (AT - 2)(iV - 3)(iV - 4)Fi£,, + (iV - 2)(iV - 3) (F^, + ISFig^,, 

+ (iV - 2) (Fi^) + SFl^ll + 18F|)^,,) + 2F^« + 2^^^)^, (26) 

where F^^^ and F^^-* are defined in Eqs. ([T8|) and (19), and F*^^-* equals 

p(3) _ rr „Pi* „P2* P3* P4 P5 P6 /oyN 

P1P2P3P4P5P6 ^ '-^^imi,£2m2/3m3/4m4/5m5/6m6 'li\mi'll2m2'lhmz'llAmA lkm^ lt-(imQ'> ' J 

where ?7^imi/2m2/3m3/4m4/5m5/6m6 cire the interaction matrix elements between symmetrized, 
normalized three-boson states. 



With the expressions (pT|-p^) it can be shown that the leading order of ^ is Further- 
more, this result remains valid for a general IBM hamiltonian with up to ra-body interactions. 
The conclusion is thus that anharmonicities can only exist in the IBM for a finite number 



of bosons. Only in the U(5) limit |^, when some of the coefficients (|2^-p6D exactly vanish, 
the spectrum can be anharmonic even for infinite N . The harmonicity of the IBM spectra 
in the limit N oo will be used as a guide to identify the number of excited bosons in a 
given state. 
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III. SCHEMATIC CALCULATIONS 



A. Dynamic symmetries 

The analysis of the three dynamical symmetry hmits of the IBM provides a good test of 
the formahsm presented in the previous section. 

In calculations carried out in the laboratory frame it is, in some cases, difficult to identify 
the character of a band 7, . . .). This becomes increasingly difficult as the number of 
phonons increases. In the intrinsic-frame calculations, on the other hand, the character of 
the band is easily identified from the structure of the state. For high the intrinsic-state 
formalism provides a good approximation to the exact calculation and it can thus be used to 
clearly identify the character of the bands in the laboratory frame. For low N the character 
of the bands in the laboratory frame can be obtained extrapolating the results of high N. 

In the following the intrinsic-state results for the three dynamical symmetry limits of 
IBM are presented. 

• U{5) limit. The hamiltonian in this limit is 

H^Sdhd^ E i^lTl-Tl. (28) 

L=l,3,4 

It corresponds to the spherical hmit (/3 — 0). There are no Goldstone bosons. The 
ground state is a condensate of s bosons (s^)^|0) and coincides with the exact (lab- 
oratory) ground state. Excited bosons correspond to S] therefore excited states are 
equal in the intrinsic and laboratory frames. A calculation in the intrinsic frame thus 
coincides with the exact one. 

• SU{2)) limit. This dynamical symmetry corresponds to the hamiltonian 

i^lTl ■ Tl. (29) 

L=l,2 

with X = i\/7/2. For prolate nuclei F^; and Ty are the only Goldstone bosons while, 
for oblate nuclei, these are F^: and F^_. The analysis of the character of the bands is 
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especially relevant in this case because the concepts of 7 and /3 bands are appropriate 
for this limit or situations close to it, i.e. well-deformed nuclei. 

The character of the bands built on one- and two-phonon excitations traditionally has 
been identified in the laboratory frame in terms of the quantum number K and the 
most symmetric representations of SU{3) [^]. In this way the f3 band corresponds to 
the representation {2N — 4, 2) K = [where (A, /i) are the Elliott quantum numbers] 
and the 7 band to {2N—A, 2) K = 2. In the case of two-phonon excitations assignments 
are already less clear. There are four possible two-phonon bands: two with K = (/?^ 
and 7^), one with K = 2 [jS'-f) and one with K = A (7^). These correspond to the 
{2N — 8,4) and the {2N — 6,0) representations of the SU{3) limit. It is not a priori 
clear, however, how the bands should be assigned to the SU{3) representations; in 
particular, this ambiguity exists for the K = bands. We now proceed to show how 
this problem can be solved using the intrinsic-state formalism. 

It is worth noting that with the definition of 7-I- and 7— bosons given in Eq. (|l^) and 
the definition of two-phonon states given in Eq. ([TTD , the states with two 7-I- or two 
7— excitations have no good K quantum number for the case of axial symmetry. In 
order to have good K in that simple linear combination correctly normalized 

can be done, 

l7ko) = ^ [l7^) + )] , (30) 

171=4) = ^ [17^) - I7' )] • (31) 

We will work with these states even in the cases in which K is considered not to be a 
good quantum number. 

In Fig. |I] we present a calculation for a particular SU (3) hamiltonian and compare 
the exact and approximate calculations for a wide range of N. In general, to identify 
the character of an excitation its energy in the intrinsic and laboratory frame should 
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be compared. This procedure is not unambiguous for low A^. Therefore, first a large 
is considered in which case the excitation energies are proportional to the number 
of excited phonons, i.e. the spectrum is harmonic (see section p.1 C| ). For ni = the 
eigenenergies of the hamiltonian (^) are 

E{L, A, /i) = -^L{L + 1) + ^(A^ + /i^ + A/i + 3A + 3/i). (32) 
o 2 

The (A,/x) contained in a symmetric representation [A^] of U{6) are 

{X,li) = {2N -An:^-6ny,2n^), n^,., = 0, 1, 2, (33) 

with 2A^ — Arix — Quy > 0. Neglecting the L{L + 1) term in (0), the excitation energies 
for the different representations are 
E{X, /x) - E{2N, 0) = -K2{Qn^ + l2ny)N + K2{Qnl + 18nJ + ISn^Uy - 3n^ - 9ny). 

(34) 

From this expression it follows that An^ = +1 corresponds to a single-phonon excita- 
tion and Auy = +1 to a double-phonon excitation. So it is clear that the one-phonon 
excitations are in the (2A^ — 4, 2) representation while the two-phonon excitations 
belong to {2N - 8, 4) or (2A^ - 6, 0). 

Finally, to fully determine the character of the states approximate and exact energies 
are compared and the structure of the state is determined in the intrinsic frame. In 
this way the mixture between the and 7|'=o excitations is found and the result 
is that (2A^ — 8,4)^ = corresponds to the combination + \I\i'k=q while 

(2A^— 6, Q)K = corresponds to ^J^P'^ — ^J^^j^^Q. If the quadrupole generator of SU (3) 
is used as E2 operator, no transitions can occur between the two bands because they 
belong to different SU{3) representations. This result is reproduced in the intrinsic 
frame since their particular composition prevents E2 transitions between them (see 
section [TVl) . 
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• 0(6) limit. The hamiltonian for this hmit is 

H = KoP^P+ J2 ^^lTl-Tl. (35) 

L=l,3 

A particular case of this hamiltonian is 

H = ■ fs, (36) 

with X = 0- This latter form is particularly useful in the study of the transition from 
SU{3) to 0(6) by varying a single parameter x- 

In this limit the concept of j3 and 7 bands is quite vague and it is therefore difficult 
to find the bandheads in the spectrum and to identify their character. To determine 
the number of excited phonons in a given state of the 0(6) limit, the harmonicity of 
the spectrum for high boson number is used (see section [II C|) . States are labeled by 
(a, z/A, T, L) and the possible values of a and r are, 

a = N -2n> 0, n = 0, 1,..., 

r = cr,cT 0. (37) 

The eigenenergies in terms of these labels are 

E{a, r, L) = Aa{a + 4) + Bt{t + 3) + CL{L + 1), (38) 

with for realistic calculations A < 0, B > and O > 0. For these parameters, 
bandheads correspond to r = and L = 0. Their energy is given by 

E{a) = Aa{a + 4) (39) 

and their excitation energy is, 

E{(y = N - In) - E{(y = N) = AA{-nN) + \A{n^ - In). (40) 

From the latter expression it follows that n represents the number of excited phonons. 
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In Fig. a calculation for a particular 0(6) hamiltonian is shown and a comparison 
of the exact and approximate results for a wide range of is presented. This limit 
corresponds to the 7- unstable vibrator and the only physical excited boson is F^. 



B. Transitional hamiltonians 

In the previous section the degree of accuracy of the method was illustrated. In this 
section schematic hamiltonians are analyzed which do not correspond to any of the three 
dynamical symmetry limits. A fixed number of bosons is used and results are plotted 
as a function of a parameter in the hamiltonian that allows to explore a wide range of 
situations, including presumably unrealistic ones. In these calculations the Goldstone bosons 
correspond to F^; and Fj, and have been removed before diagonalization. For some extreme 
values of the parameters, the ground state contains a sizeable contribution from two-phonon 
excitations. 

Figure |^a shows the results of a calculation with a quadrupole and a pairing interaction 
{H = kqP^P + K2T2 ■ T2, with X = -v^/2) for = 16. It is seen that the character of the 
different bands changes as a function of the ratio — ko/k2- Note the increase with —k,o/k,2 
of the energy of the (3 band while the 7-band energy remains approximately constant and 
the harmonic behaviour between the one-7 and two-7 bands holds throughout. This figure 



displays the same trends as Fig. 14a of Ref. p5[ which shows an identical calculation in 
the laboratory frame. Only small differences occur due to mixing with higher-order phonon 
excitations. Note that a realistic value of the ratio —k,o/h2 is around 4. 

Figure ^ shows the influence of a hexadecapole term added to a quadrupole term {H = 
K2T2 ■ T2 + ^4X4 ■ T4, with X = — /2) for A^ = 16. In the whole range of variation of the 
parameter n^j K2 there is a clear separation in energy and only small mixing between the 
one- and two-phonon states. The S\J (3) limit is recovered for K4 = 0. 

In Fig. ^ a transition from a deformed to a spherical nucleus is analyzed via the hamil- 
tonian H = Edfid + H2T2 ■ T2, with X = —V7 /2 for A^ = 16. A sharp transition occurs at a 
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particular value of the ratio —Sd/ at which point the equilibrium value of (3 drops from a 
finite value to zero. 



C. Influence of multi-phonon excitations 

Besides the mean-field approximation, another approximation made in this study is a 
truncation in the number of excited phonons that are included in the basis. This number 
is limited to two here. The previous results indicate that, in the majority of cases, excited 
states with different phonon number are decoupled. It would be of interest, however, to 
study the influence of n-phonon states with n > 2. Some idea of their influence, without 
actually extending the phonon basis, can be obtained by comparing the previous results 
(which include n = 0, 1, 2) with those where n = 0, 1. This is done in Figs. ^, ^ and |[ The 
conclusion is that the influence of two-phonon on one-phonon states is small except for some 
extreme parameter values and in regions where the two cross {e.g. for K^j ~ —4 in Fig. 



D. Simulation of triaxiality in mean field theories 

In the standard IBM, i.e. with hamiltonians up to two-body operators, a rigid triaxial 



shape cannot be obtained ||l6 ip!7| , [l9| . One way to induce triaxiality is to include a cubic 



term |]T9l into the hamiltonian. The intrinsic-state formalism supposes the presence of an 



average field, the characteristics of which should depend self-consistently on the hamiltonian. 
A simple way to investigate the consequences of triaxiality is to relax the self-consistency 
requirement and to impose by hand a triaxial mean field through a value 7 7^ 0°. It must 
be emphasized that in mean-field theories it is allowed to perturb the average field coming 
from the minimization of the energy. 

In Fig. 1^ a calculation of this type is illustrated for a simple hamiltonian varying 7. The 
value of /3 is obtained minimizing the energy for a fixed 7. For 7 > 15° the lowest band loses 
the character of ground state, i. e. this band is no longer a condensate of Fc bosons, and the 
interpretation of the levels becomes more complicated. A intriguing feature of this figure is 
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that the excitation energies of the two-phonon states are approximately constant whereas 
the excitation energies of one-phonon states decrease smoothly with 7. In other words a 
sizeable anharmonicity occurs for 7 7^ 0. For example, for 7 ~ 6° the anharmonicity q ( pT]) 
is around 0.5. This degree of anharmonicity is observed in ^^^Er and ^^^Er |I^|6|. 

This result points towards a link between anharmonicity and triaxiality. For a better un- 
derstanding of this connection, however, a fully self-consistent calculation should be carried 
out which includes three-body interactions that induce triaxiality. 

IV. ELECTROMAGNETIC TRANSITIONS 

To calculate electromagnetic transitions in the intrinsic frame one must assume that 
the total wave function can be separated into intrinsic and rotational (collective) part. 
This hypothesis is only strictly true for the SU{3) limit but it can be considered as a 
good approximation in well-deformed nuclei. The majority of deformed nuclei have axially 
symmetric shapes with an additional symmetry plane perpendicular to the nuclear symmetry 
axis. In this section we deal with this kind of nuclei. 

A general multipole operator, in the laboratory frame, of order L is related with the 
intrinsic one through 

TLMilab) = Y,V^M,fLM'{int), (41) 

M' 

where P are the rotation matrices. Laboratory and intrinsic states are also related through 
the V matrices, 

\jMK,o = 16^(^1 +^^^,) {^MK + {-iy^''vi_M^)) m). m 

where ^ [cfr. (|12D] denotes the intrinsic part of the wave function and K the projection of 
the total angular momentum on the symmetry axis. 

In our case \K^) are the intrinsic states ( p^pJ] , pOH3T| ) and reduced matrix elements in 
the laboratory frame can be related to the intrinsic matrix elements as follows: 
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J K,i\\TL{lah)\\J'K\0 = V^FTiUj'K' LK - K'\JK){Ki\fLK-K'{tnt)\K'i') 



+ {-iy'+^\j'-K'LK + K'\JK) {K^\fLK+K'\K>^>)) ] . (43) 

The last factor in r.h.s. is introduced to take into account that our intrinsic states contain 
both components ±i^. 

The reduced transition probabihty for a transition between states belonging to two ro- 
tational bands is given by 

B{E2,J'K\i' JK,0 = ^j^\{JK,^\\nilab)\\J'K',a\'- (44) 

The knowledge of the intrinsic matrix element between two rotational bands hence deter- 
mines the reduced matrix element ( ^31) from where the reduced transition probability can 
be deduced. 

Of special relevance to Nuclear Physics are E2 transitions and results given in this section 
are confined to this type of transition. The most general one-body E2 operator is defined 
in Eq. (^ (with an effective charge equal to 1). Matrix elements between the ground state 
and one-phonon states are given in Ref. [^. Likewise, matrix elements between one- and 
two-phonon can be derived. The non-vanishing ones are given by (we keep the 7-dependence 
for generality) 



(1-/3^) cos 7 - J|/3xcos27 



imoW') = \/2{N - 1) ^ ^ ' ^^^y ^ (45) 



/ ^(1 -/5^)sin7 + ^/?xsin27 



(46) 



, — sin7 + -i/|/5Ysin27 

9IT20I/57+) = Ti+F — ' ^^^^ 

; ^ COS 7 + COS 27 

|T2±2|/57+) = ^N^l ^ ^7^^ (48) 

(7+|T2o|/37+) = ^(/9|T2o|/3'), (49) 
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(7+|T2±2|/57+) = ^mf2±2\P'), (50) 



(7+|r2o|7ko) = (7+1^2017^=4) - {P\f2o\Pl+), (51) 



(7+|f2±2|7ko) = (7+|T'2±2|7k4> = (/3|f'2±2|/37+), (52) 

/? sin - + 4^ Y sin 27 - \/23 sin - + 4^ \ sin 27 
(/3^|T2±2|/3^) = (iV - 2)— + 2 , (54) 

-\/2/? sin 7 + ^ Y sin 27 
(/5^|T2o|/?7+) = \ (55) 



o - -/3cos7 + w|xcos27 
(/5^|T2±2|/37+) = ^YTF ' ^^^^ 

2/3cos7-yf/3^XCOs27 -2/3 cos 7 - ^fx cos 27 

(/37+IT20I/37+) = (iv - 2) + Y+f2 (57) 

+ Y^XCos27, 



a/2/3 sin 7 + ^/^^X sin 27 -a/2/3 sin 7 + sin 27 
1 + /92 + 1 + /32 

X sin 27 



(/37+|7^2±2|/?7+) = (A^ - 2) ^^^^^ + (58) 



V7 ' 

(/g7+|T2o|7ko) = (/37+|T2o|7l=4) = ^(/3'|T2o|/37+), (59) 

(/37+|T'2±2|7ko) = (/37+l^'2±2|7k4) = ^(/3'|T'2±2|/37+), (60) 

2/? cos 7 - a/|/9^X cos 27 /2 
(7ko|r2o|7l=o) = (7k4|7^2o|7k4) = (A^ - '^)- J^2 ^ + y 7X(cos27+ 1), 

(61) 



(7i=ol^2±2|7|=o) = (7x=4 



f2±2\ll=,) = iN-2) 



\/2/? sin 7 + y^f /^^x sin 27 



;Xsin27 



(62) 



In order to see the general behaviour of the calculated transitions, in Fig. ^ we present a 
calculation for the B{E2)^s as a function of x fo^' the hamiltonian H = K2 {T2 ■ ^2 "~ P^'P^- 
In (a) the calculation is done directly in the lab system while in (b) the intrinsic wave 
functions and Eqs. (^3| - ^ ) are used. The main characteristics of the calculated B{E2)^s are 



similar to those obtained in the lab in Ref. p5| although in that case only ground state and 
one phonon excitations are included. 



V. CONCLUSIONS 

In this paper a formalism is presented for performing calculations with up to two phonon- 
excitations in the intrinsic frame. This formalism yields good agreement with exact results 
and allows an easy interpretation of them. The calculations can be easily extended to include 
higher multi-phonon excitations and/or higher-order interactions. 

From these results it is deduced that for a large number of bosons the IBM is a 
harmonic model even if n-body interactions are included. A sizeable degree of anharmonicity 
can be obtained only for a finite number of bosons combined with a three-body term in the 
hamiltonian. This result points towards a link between triaxiality and anharmonicity in the 
IBM. 
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FIGURES 

FIG. 1. Energy levels as a function of the number of bosons N for the SU{3) limit. The 
hamiltonian H = k2T2 ■ T2 is used. Full lines correspond to the intrinsic-frame calculations and 
dashed lines to the laboratory-frame calculations. 

FIG. 2. Energy levels as a function of the number of bosons N for the 0(6) limit. The hamil- 
tonian H = koP^P is used. Full lines correspond to the intrinsic- frame calculations and dashed 
lines to the laboratory-frame calculations. 

FIG. 3. Energy levels and composition of states as a function of the ratio —hq/k2 for the 
hamiltonian H = kqP'^P + K2T2 ■ T2, with x = — \/7/2 and = 16 for up to (a) two phonons and 
(b) one phonon. 

FIG. 4. Energy levels and composition of the states as a function of the ratio —K4/K2 for the 
hamiltonian H = K.2T2 • T2 + • T4, with x = —V^ /2 and = 16 for up to (a) two phonons 
and (b) one phonon. 

FIG. 5. Energy levels and composition of the states as a function of the ratio —£d/K-2 for the 
hamiltonian H = Edfid + i^2T2 ■ T2, with x = — \/7/2 and = 16 for up to (a) two phonons and 
(b) one phonon. 

FIG. 6. Excitation energies and composition of the states as a function of 7 in a 7-perturbed 
mean-field calculation for the hamiltonian H = K2{T2 • T2 — P'^P), with x = ~V^/2 and N = 16. 

FIG. 7. B{E2) values for selected transitions between ground, one- and two-phonon excited 
states for a hamiltonian H = K2 (T2 ■ T2 — P^P^ as a function of x for a system with = 16. In 
(a) the calculation is done directly in the lab system and in (b) the intrinsic wave functions and 
Eqs. (||-5D are used. 
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